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CROSS-REFERENCE TO RELATED APPLICATIONS 

[0001] This patent application claims priority from Provisional Patent Application No. 
60/415,510 filed October 1, 2002 for "METHOD FOR RAPID SCATTERING 
SIMULATION OF OBJECTS IN. IMAGING VIA DOMAIN DECOMPOSITION AND 
SPECTRAL MATCHING," and is related to co-pending application serial no. 10/241,242 
filed September 10, 2002 for "CHARACTERIZING ABERRATIONS IN AN IMAGING 
LENS AND APPLICATIONS TO VISUAL TESTING AND INTEGRATED CIRCUIT 
MASK ANALYSIS," which are incorporated herein by reference for all purposes. 

BACKGROUND OF THE INVENTION 

[00021 This invention relates generally to optical imaging systems, and more particularly 
the invention relates to simulation and modeling of electromagnetic scattering of light in such 
imaging systems. The invention has applicability to optical proximity correction in photo 
masks and mask image inspection, but the invention is not limited thereto. 
[0003] The general components of an optical lithography tool, shown schematically in the 
diagram of Figure 1 are the illumination system, the projection system, the photomask (also 
called reticle), and the photoresist, spun on top of a semiconductor wafer. The operation 
principle of the system is based on the ability of the resist to record an image of the pattern to 
be printed. The mask, already carrying this pattern, is flooded with Hght and the projector 
forms an image of all mask patterns simultaneously onto (and into) the resist. The inherent 
parallelism of this process is the main reason why optical photolithography is favored over 
any other lithography, since it facilitates a very high throughput of 30-120 wafers per hour. 
The Ught intensity distribution on top of the resist surface is commonly referred to as aerial 
image. The resist itself is a photosensitive material whose chemical composition changes 
during light exposure. The pattern is thereby stored in form of a latent (bulk) image within 
the resist. After exposure has occurred, the resist is developed by means of a chemical 
process that resembles the process of developing photographic film. After development, the 
exposed parts of the resist remain or dissolve depending on its polarity (negative or positive, 
respectively). The end-result of the lithography process is a more or less exact (scaled or not) 
replica of the mask pattern on the wafer surface that will play the role of a local protective 



UCALP016 



Page 1 



layer (mask) for subsequent processing steps (etching, deposition, implantation). 
[00041 The role of the illumination system is to deliver a light beam that uniformly trans- 
illuminates the entire reticle. See Figure 2(a). It typically consists of various optical 
elements, such as lenses, apertures, filters and mirrors. The Ught source is responsible for 
generating very powerful and monochromatic radiation. Power is necessary because it is 
directly related to throughput. Monochromaticity is important because high quality refractive 
(or reflective, in the case of EUV lithography) optics can only be fabricated for a very narrow 
illumination bandwidth. State-of-the-art optical lithography tools employ excimer lasers as 
their light source. Deep Ultra Violet (DUV) lithography is the term used for lithography 
systems with illumination wavelengths ?.=248nm (excimer laser with KrF), X=193nm 
(excimer laser with ArF) and ^=157nm (excimer laser with F2). The successful development 
of current and future optical photolithography technologies is hinged upon research advances 
in both excimer laser technology and novel materials that possess the required properties 
(high optical transmission at DUV wavelengths, thermal properties, stability after heavy 
DUV radiation exposure) by which the optical elements of the system will be made. 
[00051 All illumination systems in optical projection printing tools are designed to 
provide what is known as KOhler illumination. By placing the source or an image of the 
source in the front focal plane of the condenser column, the rays originating from each source 
point illuminate the mask as a parallel beam, as seen in Figure 2. Here, each source point 
emits a spherical wave that is converted by the illumination system into a plane wave incident 
on the object (photomask). The angle of incidence of the plane wave depends on the location 
of the source point (a,/3) with respect to the optical axis (0,0). Each parallel beam is a plane 
wave whose direction of propagation depends on the relative position of the source point with 
respect to the optical axis. Nonuniformity in the brightness of the source points is averaged 
out so that every location on the reticle receives the same amount of illumination energy. As 
we will see in subsequent Sections Kohler illumination can be modeled in a concise 
mathematical way. 

[00061 In addition to dose uniformity, the lithography process should also maintain 
directional uniformity such that the same features are replicated identically regardless of their 
orientations. The shape of the Ught source is therefore circular (or rotationally symmetric) in 
traditional optical lithography, although this is not true for certain advanced illumination 
schemes such as quadrupole illumination, where directional uniformity is sacrificed in order 
to maximize the resolution of features with certain orientations. 
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[0007] The coherence of the light source is another important attribute. Temporal 
coherence is usually not a big concern, since the narrow bandwidth of excimer lasers implies 
high temporal coherence. Spatial coherence (or just plain coherence) on the other hand is 
always carefully engineered and, in most cases, adjustable. Using special scrambling 
techniques, the Ught emitted from any point of the source is made completely uncorrelated 
(incoherent) to the light emitted from every other point. However, Ught gathers coherence as 
it propagates away from its source. The frequently quoted partial coherence factor a is a 
characteristic of the illumination system and is a measure of the physical extent and shape of 
the light source. The larger the light source, the larger the partial coherence factor, and the 
Ught source has a lower degree of coherence. In the limit of an infinite source, imaging is 
incoherent and a - oo . On the other hand, the smaller the light source, the smaller the partial 
coherence factor, and the higher the degree of coherence. Unaging with a point source is 
fully coherent and a = 0. Note that a point source in a Kohler illumination will result in a 
single plane wave illuminating the mask and the angle of incidence of this wave depends on 
the relative position of the point source witii respect to the optical axis. For partial coherence 
factors between zero and infinity, imaging is partially coherent. Typical partial coherence 
factors in optical litiiography range from 0.3 to 0.9. 

[0008] The projection system typically consists of a multi-element lens column (up to 30- 
40 lenses) that may also have apertures, filters, or other optical elements, and it is a marvel of 
engineering precision in order to be able to reliably project images with minimum dimensions 
on the order of lOOnm for state-of-the-art systems. One of the main reasons for the required 
high precision is control of the aberrations or deviations of the wavefront from its ideal shape. 
Two relevant parameters of the projection system are the numerical aperture, NA, and the 
reduction factor, R. The numerical aperture is, by definition, the sine of the half-angle of the 
acceptance cone of light-rays as seen from the image side of the system. The ratio of image- 
height to object-height is, by definition, tiie magnification factor M of the system. The 
inverse of the magnification factor is the reduction factor, R. Since a typical system in 
photolithography projects at the image plane a scaled-down version of tiie object (mask), M 
is less tiian 1 and R is greater than 1. State-of-the-art systems currently have reduction 
factors of R=4 or 5. Note that two numerical apertiires exist in the projection system, namely 
NAi (or simply NA) and NAo, which refer to the half-angle of the acceptance cone as seen 
from the image side and from the object (mask) side, respectively. They are related through 
the reduction factor as follows: 
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Equation 2-1. NA 



[00091 For a circularly shaped light source, the partial coherence factor a mentioned 
above is related to the numerical apertures of both the projection system and the illumination 
system. Specifically, ct is given by: 



_ NA 

Equation 2-2. ^ ~ NA 



p 



where NAc is the numerical aperture of the condenser lens (illumination system) and NAp is 
the numerical aperture of the projector lens. Some confusion arises firom the fact that, in the 
above equation, the reduction factor of the imaging system is implicitly taken into account. 
Figure 2(b) clarifies the situation by showing simplified diagrams of two optical systems with 
parameters NA=0.5, a=0.5 and R=5 or R=l. 

[00101 The photomask, also called reticle, carries the pattern to be printed at a given 
Uthography processing step. The masks of integrated circuits having large die-sizes or 
footprints (that is, occupying large areas on the semiconductor wafer), typically carry just one 
copy of the chip pattern. A matrix of several chip patterns is contained in one mask 
whenever the chip size permits. Note that the mask is drawn R tunes the actual size on the 
semiconductor wafer, since the dimensions of the circuit will be scaled down by the reduction 
factor, R. For this reason it is not sufficient to just provide feature sizes, since it may not be 
immediately obvious firom the context whether these are photomask (object) or resist (image) 
sizes. A typical convention for distinguishing photomask feature sizes fi-om resist feature 
sizes is to include in parenthesis the reduction factor R. For example, a 600mn (4X) line has 
a size of 600nm on the mask, and would produce a 600mn/4=150nm line if used in a 4X 
imaging system. Similarly, a 130mn (IX) line refers to the size of a line at the image (wafer) 
plane and would result from the printing of a 130nm line on the mask for a system with R=l, 
or a 520nm (4xl30nm) line on the mask for a system with R=4, or a 1 .3^im (10xl30nm) line 
on the mask for a system with R=10. 

[0011] Depending on their operation principle, photomasks can be divided into two broad 
categories: conventional binary or chrome-on-glass (COG) masks and advanced phase- 
shifting masks (PSM). 

[0012] A binary or COG mask consists of a transparent substrate (mask blank), covered 
with a thin opaque film that bears the desired pattern. Light can either pass unobstructed 
through an area not covered by the opaque film or be completely blocked if it is incident on 
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an area that is protected by the film. This binary behavior of the transmission characteristic 
of the mask is responsible for its name. The mask blank for DUV lithography typically 
consists of fused silica glass that has excellent transmission at X=248nm and somewhat 
poorer but acceptable transmission at X,=193mn and X=157mn. The opaque fihn is typically 
on the order of lOOnm thick and has a chromium (Cr) composition. 

[0013] Adding phase modulation to the photomask can profoundly increase the attainable 
resolution. This is the principle followed by phase-shifting masks (PSM), which employ 
discrete transmission and discrete phase modulation. There are many different flavors of 
PSMs depending on the way that the phase modulation is achieved. One of the most 
promising PSM technologies is what is known as alternating phase-shifting mask (alt. PSM, 
or APSM). Here are cut-planes of geometry of a binary (COG) mask (a) and an alternating 
phase-shift mask (b). The ideal electric field distribution for the binary mask (c) leads to a 
poor image intensity distribution (e) at the image plane, whereas the ideal electric field 
distribution for the alt. PSM (d), because of destructive interference, leads to a robust image. 
The principle of an alt. PSM is compared with that of a binary mask in Figure 3. The center 
line is bordered by transmitting regions with 1 80° phase difference on an alt. PSM and by 
clear areas of the same phase on a binary mask. The phase difference on the alt. PSM leads 
to destructive interference, resulting in a sharp dark image. The binary mask image is not as 
sharp because ofthe lack of phase interaction. The 180° phase difference is created by 
etching trenches, also called phase-wells, into the fiised silica substrate during the alt. PSM 
fabrication process, which is now more complex than the COG fabrication process. The 
difference in the amount of material removed detch is such that the path length difference 
between Ught passing through the different phase regions is half ofthe wavelength in air. 
[00141 Sub-wavelength lithography, where the size of printed features is smaller than the 
exposure wavelength, places a tremendous burden on the lithographic process. Distortions of 
the intended images inevitably arise, primarily because ofthe nonlinearities ofthe imaging 
process and the nonlinear response ofthe photoresist. Two ofthe most prominent types of 
distortions are the wide variation in the linewidths of identically drawn features in dense and 
isolated enviromnents (dense-iso bias) and the line-end puUback or line-end shortening (LES) 
from drawn positions. The former type of distortion can cause variations in circuit timing 
and yield, whereas the latter can lead to poor current tolerances and higher probabilities of 
electrical failure. 

[0015] Optical proximity correction or optical proximity compensation (OPC) is the 
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technology used to compensate for these types of distortions. OPC is loosely defined as the 
procedure of compensating (pre-distorting) the mask layout of the critical IC layers for the 
lithographic process distortions to follow. This is done with specialized OPC software, hi 
the heart of the OPC software is a mathematical description of the process distortions. This 
description can either be in the form of simple shape manipulation rules, in which case the 
OPC is referred to as "rule-based OPC," or a more detailed and intricate process model for a 
"model-based OPC." The OPC software automatically changes the mask layout by moving 
segments of line edges and adding extra features that (pre-) compensate the layout for the 
distortions to come. Although after OPC has been performed the mask layout may be quite 
different than the original (before OPC) mask, the net result of this procedure is a printed 
pattern on the wafer that is closest to the IC designer's original intent. 
[00161 In the early 1990's the problem of OPC for large mask layouts was addressed 
formally as an optimization problem. The size of such a problem is formidable, but through 
introduction of appropriate constraints - induced primarily from mask fabrication 
constraints, local optimal mask "points" were successfully demonstrated. This research has 
led to commercially available software tools that perform OPC on a full-chip scale. All three 
approaches rely heavily on speedy calculations of the image intensity at selected points of the 
image field. Although the methods by which they achieve this appear seemingly different, 
they are nevertheless the same, in the sense that they all rely on a decomposition of the kernel 
of the imaging equation for partially coherent light. 

[0017] Domain decomposition techniques, where a large electromagnetic problem is 
broken up into smaller pieces and the final solution is arrived at by synthesizing (field- 
stitching) the elemental solutions, have been proposed for the study of one-dimensional 
binary (phase only) diffractive optical elements. Others working on the same problem have 
demonstrated how to create and use a perturbation model for binary edge-transitions based on 
the product of the ideal, sharp transition and the continuous field variations in the vicinity of 
the edge. Independently, a similar technique to field-stitching has been proposed for the 
simulation of two-dimensional layouts of advanced photomasks (altemating PSM, masks 
with OPC) that properly models interactions &om neighboring apertures and fiirthermore 
takes advantage of the spectral properties of the diffracted fields to come up with a compact 
model for the edge-transitions. 

[0018] The present invention is directed to an improved methodology for OPC on optical 
masks. 



UCALP016 



Page 6 



SUMMARY OF THE INVENTION 

[00191 In accordance with the invention, a complex two-dimensional layout of a 
photomask or other three-dunensional object is systematically decomposed into a finite 
number of elementary two-dimensional objects with the abiUty to cause one-dimensional 
changes in light transmission properties. An algorithmic implementation of this can take the 
form of creation of a look-up table that stores all the scattering information of all two- 
dimensional objects needed for the synthesis of the electromagnetic scattered field Scorn the 
original three-dimensional object. 

[0020] More particularly, the domain is decomposed into edges, where pre-calculated 
electromagnetic field from the diffraction of isolated edges is recycled in the synthesis of the 
near diffracted field from arbitrary two-dimensional diffracting geometries. In this invention, 
the sum of the complex fields is taken rather than the product due to the linearity of the 
Kirchoff-Fresnel diffraction integral, which makes the decomposition method possible. 
[0021] The invention has particular applicability in die-to-database inspection where an 
actual image of a mask is compared with a synthesized image that takes imaging artifacts of 
comers, edges and proximity into account. A very important current challenge in 
synthesizing images is the correct evaluation of elecfromagnetic effects at edges of phase- 
shifted mask features. These electromagnetic eff^ects are a sfrong fimction of the edge 
topography of the nonplanar features infroduced during the mask writing process. The 
inspection of every featiire on a 50M transistor chip must be accomplished in an inspection 
time on the order of 5 hours for feature sizes of 1 OOnm. Futiire inspection rates are expected 
to hit 200M inspection pixels per hour, which would allow a mask with 50nm features and 
200M transistors to be inspected at a pixel size of 25nm in about 5 hours. 
[00221 Another application is optical proximity correction which consists of evaluating 
the image of every feature on a mask and improving it by introducing edge shifts and 
iteratively adjusting the amounts of these shifts. A typical OPC correction process on a 50M 
transistor mask takes about a full day. Since the number of iterations in the OPC shifts is on 
the order of 5 the M\ mask image is computed on average roughly at the same speed of the 
die-to-database inspection, or 50M transistors every 5 hours for the 100 nm generation and 
200M transistors per hour for the 50 nm generation. 

[00231 The invention and objects, features and appUcations thereof will be more apparent 
from the following detailed description and appended claims when taken with the drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[00241 Fig. 1 is a general diagram of an optical photolithography tool. 

[00251 Fig. 2(a) illustrates the illumination system of Fig. 1 with Kohler's method. 

[00261 Fig. 2(b) illustrates numerical apertures of illumination and projection lenses in 

the tool of Fig. 1. 

[00271 Figs. 3(a)-3(f) are a comparison of operational principles of an alternating phase 

shift mask (alt. PSM) and a binary (COG) mask as used in the invention. 

[00281 Fig. 4 illustrates light decomposition of an opening into two edges. 

[0029] Figs. 5(a)-5(e) illustrate an edge domain decomposition method in accordance 

with the invention. 

[0030] Figs. 6(a)-6(e) illustrate an edge domain decomposition method applied on an 
isolated line. 

[0031] Fig. 7 illustrates an edge domain decomposition method applied on an arbitrary 
polygon. 

[00321 Figs. 8(a)-8(d) illustrate dependence of edge scattering on edge profile and 
polarization. 

[0033] Figs. 9(a)-9(d) illustrate complex differences of true edge-scattering firom ideal 
and sharp transition. 

[00341 Figs. 10(a)-10(d) illustrate the spectra of complex differences shown in Figs. 9A- 
9D. 

[0035] Fig. 1 1 illustrates edge domain decomposition on a large mask layout. 
[0036] Figs. 12(a)-12(c) illustrate edge domain decomposition applied on an isolated 

square hole (scattered field). 

[0037] Figs. 13(a)-13(c) illustrate edge domain decomposition applied on an isolated 
square hole (spectrum). 

[00381 Figs. 14(a)-14(c) illustrate edge domain decomposition applied on an isolated 
square island (scattered field). 

[00391 Figs. 15(a)-15(c) illustrate edge domain decomposition applied on an isolated 
square island (spectrum). 

[0040] Fig. 16 illustrates edge domain decomposition on a large mask layout (images). 
[0041] Figs. 17(a), 17(b) illustrate matched bandwidth approximation to edge diffi-action. 
[0042] Figs. 18(a), 18(b) illustrate matched bandwidth approximation applied to the hole 
of Fig. 12. 

[0043] Figs. 1 9(a), 19(b) illustrate defect assessment using prior art techniques. 
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[00441 Fig. 20 illustrates domain decomposition in accordance with an embodiment of 
the invention when a defect is present. 

[00451 Figs. 21(a), 21(b) illustrate layout and amplitude of scattered field below the mask 
at the observation plane. 

[00461 Figs. 22(a)-22(d) illustrate a 90 degree phase defect present in the center of the 
bottom-left hole of a mask. 

[00471 Figs. 23(a)-23(d) illustrate a defect present in the center of the bottom-right hole 
of a mask. 

[00481 Figs. 24(a)-24(b) are a comparison of the decomposition method with the rigorous 

determined scattering. 

[00491 Figs. 25(a)-25(d) are plan views of defective masks of Figs. 22 and 23. 
[00501 Figs. 26(a)-26(b) are spectra plots (magnitude) of errors of Figs. 22 and 23. 

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS 
[00511 Adam and Neureuther, "SIMPLIFIED MODELS FOR EDGE TRANSITIONS IN 
RIGOROUS MASK MODELING," 26* Annual Intemational Symposium on 
Microlithography, March 1, 2001, discloses domain decomposition methods (DDM) that 
allow the division of a larger mask diffraction problem into a set of constituent parts or single 
openings in the mask. 

[00521 The present invention provides an extension of the disclosed domain 
decomposition methods (DDM). A domain decomposition method based on edges uses pre- 
calculated electromagnetic field from the diffraction of isolated edges which is recycled m 
the synthesis of the near diffracted field from arbitrary two-dimensional diffracting 
geometries. Although at first glance it might be tempting to discard edge-decomposition 
techniques as inaccurate for masks with large vertical topography and features with small 
■ lateral dimensions, it will be shown that for a large set of practical situations, they are 
accurate. A key difference in this method from the prior art is that the sum of the complex 
fields is taken instead of the product. The reason for this is in the linearity of the Kirchhoff- 
Fresnel diffraction integral, and this linearity makes possible the decomposition method. 
[00531 A natural extension of the DDM that is also based on the linearity of the 
Kirchhoff-Fresnel diffraction integral is shown in Figure 4. An opening with transmission T 
and phase d, in an infinitely thin, opaque screen can be decomposed by way of the linearity of 
the Kirchhoff-Fresnel diffraction integral into two edges. The uniform plane wave 
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illumination is subtracted to restore the light level everywhere. The ID mask layout is first 
decomposed into single opening masks, as noted above, and subsequently each opening is 
decomposed into two edges and a uniform illumination field as in Figure 4. Such a 
decomposition can be proven to be exact for infinitely thin and perfectly conducting screens 
through the use of the electromagnetic form of Babinet's principle, but it should break down 
if applied to photomasks with "thick" vertical structures (compared to X) and not perfectly 
absorbing screen materials. Also, because of the reduced abiUty of Ught to penetrate through 
small (compared to X) openings, it is expected to break down for small mask features. This 
decomposition method will be hereafter referred to as the edge domain decomposition 
method or edge-DDM. The application of the edge-DDM to a single isolated space and an 
isolated line are described next, followed by a systematic evaluation of the limits of 
applicability of the method. 

[0054] Consider the isolated mask space (opening) shown in Figures 5(a)-5(f). The 
amplitude of the complex Ey-field under Ey (TE) illumination is shown in Figure 5(c) and 
the scattered field across the observation plane is shown in Figure 5E. According to Figure 4, 
the field across the observation plane can be alternatively obtained by using the scattered 
field from an edge that has the same vertical profile as the space, as shown in Figures 5(b) 
and 5(d) The resulting field from this edge-DDM is overlaid on the plot of Figure 5E along 
with the field from the edge-scattering. The normalized mean square errori (NMSE) between 
the r-mask and the edge-DDM is 0.19%. In Figure 5(f) the spectra of the r-mask and edge- 
DDM are compared and it is seen that they differ by O.I30/0 in a NMSE sense if the whole 
spectrum of propagating waves is considered and by merely O.OO70/0 within the collection 
ability of a projection system with NA=0.8 and R=4. Since only the portion of the spectrum 
that is collected by the projection system contributes to the image formation, the exceptional 
accuracy of the edge-DDM is this regime is more than adequate for accurate image 
simulations. 

100551 Similar steps can be followed for the isolated line shown in Figure 6(a) 
illuminated with Ex (TM) polarized light. Note that it is not typical to have a line as shown 
here, where in both sides the glass has been etched. Normally only one side of the line is 
phaseshifted. but this example simply aims to introduce the principle of the edge-DDM. A 
line with one side phase-shifted and the other left intact can still be decomposed via the edge- 
DDM using the diffraction fields from the two different edges. The amplitude of the complex 
Ex field everywhere in the domain of Figure 6(a) is shown in Figure 6(c) and for the edge of 



UCALP016 



Page 10 



Figure 6(b) in Figure 6(d) respectively. The true (r-mask) scattered field is compared with 
the synthesized field with edge-DDM in Figure 6(e) and the spectra are compared in Figure 
6(f). The error (NMSE) is seen to be small in both the near field and spectra plots and 
specifically the error of the through-the-lens (TTL) spectrum is just 1.5E-5. 
[00561 In the examples of Figure 5 and Figure 6 the edge-DDM is seen to be accurate 
compared with the simulation of the exact mask structure. Although the vertical mask 
topography (180° « 170nm « IX) is comparable to the wavelength, the fact that the lateral 
dimensions of the space and the line are relatively large compared with the wavelength (CD = 
400mn > 2 x 193mn = 2 x X) renders the edge-DDM accurate. One expects that as the lateral 
sizes become smaller and the vertical mask topography even larger the edge-DDM will 
eventually break down. The limits of the edge-DDM were investigated through exhaustive 
shnulation of isolated spaces and lines on the mask with various etched depths and amounts 
of undercut. Hie sidewall angles were always kept vertical and no rounding at the bottom of 
the etched wells was included. The mask CD was varied in the set {2^m, 1 ^m, 0.6nm, 
0 4^ml, 0.2^m(~a), O.^m}, the etched depth in {0», 90», 180», 270», 360'>} and the amount 
of underetch in {Omn, 25mn, 50mn, 75mn, lOOmn}. All 2 (spaceAine) by 2 (TE/TM 
illumination) by 6 (CD sizes) by 5 (etched depths) by 5 (underetch) = 600 cases were run and 
the rigorous (r-mask) scattered field was automatically compared with that obtained firom 
edge-DDM. The results of the investigation were spectacular. The edge-DDM can readily 
achieve better than 1% accuracy (in a NMSE sense) in the near field for CD=400mn (4X) or 
larger regardless of depth, underetch and polarization. The accuracy in the TTL-spectrum is 
sufficient in those cases even for inspection (R=l) simulations. For lines, the edge-DDM 
maintains excellent accuracy in the TTL-spectrum with R=4 down to a mask CD of lOOmn 
(~0 5X) with 90- of etched depth and Omn of undercut, or mask CD of 200nm (~U) with 270- 
of etched depth and 25mn of undercut. For spaces, the accuracy of the method is acceptable 
at least down to mask CD of 200mn (-IX) with 270- of etched depth and 50mn of underetch. 
[0057] Current state-of-the-art photomasks and processes utilize sub-resolution assist 
features that are never less than one wavelength in size on the mask. This clearly implies that 
the edge-DDM has the potential to accurately simulate the most advanced photomask 
technologies (OPC, alt. PSM, OPC on alt. PSM). Again, the true power of the method lies on 
its application in 2D mask layouts. However, one key difference with the DDM of mask 
openings is the following: only one rigorous 2D simulation of an edge suffices for the 
scattered field reconstruction of any arbitrary size opening on the mask. The DDM of 
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openings requires a separate simulation of every different mask size. Next, the application of 
the edge-DDM in 2D layouts is developed. 

[00581 The application of the edge-DDM in 2D layouts is straightforward. A mask 
layout comprises of a large number of edges positioned at different locations and having 
different orientations. No matter how complicated the mask technology, there are usually 
only a small number of different types of edges present in the layout. For example, a single 
exposure 0°/90»/270° alt. PSM has five types of edges: i) Cr-layer/O" edge, ii) Cr-layer/90° 
edge, iii) Cr-layer/270» edge, iv) 0° /90° edge and /270' edge. Depending on the 
orientation of each edge in the layout it "sees" and responds to the incident field differently. 
If the illuminating field is a TE (Ey) normally incident plane wave and the edge is oriented 
along the y-axis it "sees" TE illumination, but if it is oriented along the x-axis it "sees" TM 
illumination. This is illustrated in Figure 7. Depending upon the orientation of each edge in 
the layout it "sees" and responds to the incident field differently. Here the illuminating field 
is a TE (Ey) normally incident plane wave and the edges that are oriented along the y-axis 
"see" TE illumination (parallel to the edge), whereas the edges that are oriented along the x- 
axis "see" TM illumination (perpendicular to the edge). A separate 2D edge simulation is 
required to capture the different response of the edge to different field polarizations. The 
response of each edge to its respective illumination is taken into account rigorously firom a 
pre-stored 2D edge-diffraction simulation. However, comer effects are effectively ignored, 
since the finite extent of each edge is not rigorously taken mto account. Instead, the scattered 
field at the end points of every edge is abruptly terminated (truncated) in a perfect square- 
wave fashion to the field value of the kmask model. Comer effects in typical imaging 
situations are mapped at the extremities of the spectra and do not contribute to the image 
formation. Owing to this observation, the edge-DDM applied in 2D layouts is successM m 
accurately capturing the true electromagnetic behavior of edges in a rapid mamier, while the 
insignificant comer effects are safely ignored. 

100591 Slightly different embodiments of the edge-DDM applied in 2D layouts can be 
devised. One that is powerfiil and general in its appUcation is described next. Then, 
scattering results from various edge profiles that would be encountered in alt. PSMs are 
presented for completeness. Next, application of the edge-DDM in simple 2D mask layouts 
is presented and the results of the method are compared with fiiUy rigorous r-mask models m 
order to establish the validity of this approach. Finally, a large portion of a real layout from a 
single exposure 07907270° alt. PSM that is too big for 3D rigorous simulation is simulated 
via the edge-DDM and the aerial image is compared with the simple k-mask model. 
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[0060] First, every distinct edge that is present on the layout is pre-simulated with all 
possible illumination directions (field polarizations) that are required, based on the 
orientations of the edge encountered within the layout. For simplicity only Manhattan-type 
layouts are considered, where the edges are oriented along either the x- or the y-axis and 
require electromagnetic simulation under both parallel and perpendicular to the edge 
polarizations. The extension to non-Manhattan layouts although much more computationally 
laborious is straightforward. Subsequently, each distinct mask layer is broken up into a set of 
mutually disjoint (non-overlapping) rectangles. For example, a plain binary mask consists of 
two layers, namely the absorption layer (typically Cr-based) and the clear (uncovered) layer, 
whereas a more advanced 0«'/90°/270» three-phase alt. PSM consists of four layers, namely the 
absorption layer and the 0'' clear, 90" clear and 270° clear layers. Next, looping through all 
rectangles of each layer, the type of each one of the four edges of the rectangle is determined 
based on what its neighboring layer is. Trivial edges (when the neighboring rectangle is of 
the same type) are discarded. It is also possible that only part of a single edge from a 
rectangle is neighboring with one layer and the rest is neighboring with a different layer (or 
different layers). In that case the edge is broken up in multiple edges, in a way that tracks the 
neighboring layer. What is left at the end of this process is a set of all non-trivial edges that 
are present in the layout (location, size and orientation) and their type. i.e. which are the two 
mask layers on each side of the edge. Given the mask layout, the k-mask (Kirchhoff-mask) 
model is readily available. Finally, the respective (complex) difference of the true edge- 
scattering from the ideal and sharp k-mask model is added to all non-trivial edges 
accordingly. This revised mask model that results from the edge-DDM includes accurate 
information about the electromagnetic scattering from the edges and it is a quasi-rigorous 
mask model (qr-mask) in the same sense that the term was used where only 2D mask 
simulations were used to approximate a 3D electromagnetic problem. 

[0061] The above algorithm was implemented in a MATLAB environment. The part that 
deals with the polygon and edge extraction can probably be judged as rudimentary and most 
likely cannot compete in speed with more sophisticated implementations in production 
caliber CAD software, but it suffices for the purposes of this work, that is, to demonstrate the 
proof of concept of the edge-DDM and the feasibility of including accurate edgescattering 
information in rapid aerial image simulations. 

[0062] Examining the departure of the scattered field from an edge from its ideal, sharp 
edge-transition and how it depends on the profile of the edge and the incident field 
polarization can provide intuition of the physical mechanisms involved. 



UCALP016 



Page 13 



[0063] Figure 8 shows the amplitude of the (near) scattered field from an edge with 
various profiles. The ideal, sharp edge-transition is also shown on every graph for 
comparison. In Figure 8(a) and (b) the edge-scattering for different etched depths (0°, 180° 
and 360°) and vertical sidewalls with no undercut is shown for TE (Ey) and TM (Ex) 
illumination respectively. The effect of undercut is depicted in Figure 8(c) and (d), where the 
edgescattering for different amounts of undercut (Omn, 25nm, 50mn and 75mn) and vertical 
sidewalls with 1 80» etched depth is shown for TE and TM illumination respectively. It is 
interesting to note that increasing the etched depth is "making the edge response slower", 
meaning that the near scattered field from the edge requires more distance to reach its "on" 
value from its "off value, although this effect is more pronounced for the first increment 
from 0° to 180° rather than from 180° to 360». Similarly, increasing the amount of undercut 
appears to "make the edge response faster" and furthermore this behavior appears to be 
linear, since for every increment of the undercut the edge-transition gains a constant distance 
in its "off-on" behavior. 

[00641 It is instructive to isolate the difference of the true edge- scattering from the ideal, 
sharp edge-transition. This is done in Figure 9, which depicts the amplitude of this (complex) 
difference for all respective plots of Figure 8. Here the amplitude of the complex difference 
of the true electromagnetic field scattered from an edge from the ideal, sharp transition (step- 
function) is shown in Figures 9(a)-(d) for all edge profiles and illumination polarizations of 
Figure 8 respectively. Figures 9(a) and (c) are for TE polarization and Figures 9(b) and (d) 
are for TM polarization. The spatial extent of the difference grows with increasing etched 
depths, as seen in Figures 9 (a) and (b), and shrinks with increasing undercuts, as seen in 
Figure 9(d) for TM polarization. The behavior in Figure 9(c) appears somewhat erratic, but if 
the phase is also taken into account, it is shown in Figure 10 that this is also the case for TE 
polarization. Observe that, indeed, by increasing the etched depth the difference has a larger 
spatial extent, whereas increasing the amount of undercut confmes the spatial extent of the 
difference. This is crucial, because it is exactly the spatial extent of the difference that affects 
its spectral distribution. If the difference is well confined, its spectral disfribution is spread 
out to higher spatial frequencies and the lower frequency content is small. In such a case the 
complex difference can be neglected and the approximation of the ideal, sharp transition is 
adequate. 

[0065] However, when the difference has a larger spatial extent, then its lower frequency 
content becomes significant and it camiot be neglected. iWs is shown in Figure 10. which 
depicts the magnitude of the spectra of all respective difference disfributions of Figure 9. 
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Here the spectra of the complex differences shown in Figure 9 are depicted here respectively. 
Figures 10(a) and (c) are for TE polarization and Figures 10(b) and (d) are for TM 
polarization. Observe that the through-the-lens (TTL) parts of these spectra, which for ^ 
typical lithographic imaging situations (R=4, NA~0.7-0.8) extent up to ~ ± 0.008 (mn) ' are 
smaller for smaller etched depths and larger undercuts, regardless of polarization, as expected 
from the observations on the plots of Figure 8 and Figure 9. Observe that the TTL spectrum 
of the difference attains larger values when the etched depth increases and smaller values 
when the amount of undercut increases, as expected from Figure 8 and Figure 9. Remember 
that the normalization used in spectra plots is that a constant (DC) value of 1 (clear field) 
would produce a Dirac 6-function at zero frequency with amplitude of 19.4 

(V^ = V>/^"19-4")- 

[0066] The application of the edge-DDM on an isolated square hole and an isolated 
square island (post) and comparison of the results with fixUy rigorous 3D mask simulations 
(rmasks) are presented here. The edge profile has a glass etched depth of 180» and an 
undercut of 50mn. The dimensions of both the hole and the island are 400mn x 400mn (4X), 
corresponding to approximately 2x2 wavelengths, for X,=193nm. 
[00671 Figure 12(a) and (b) depict the amplitude of the scattered field across the 
observation plane below the hole obtained via the edge-DDM and via M\y rigorous 3D 
simulation, respectively. The amplitude of their difference (error of the edge-DDM) is shown 
in Figure 12(c). Locally the amplitude of the error is seen to reach a discouragmg level of 
almost 50% of the clear field value. The normalized mean square error is however more 
contained and is only about 4.9% if calculated within a 3^lm by 3^m area or 3.6% if 
calculated within the 0.4^m by 0.4^m open area of the hole, as shown in Figure 12. From 
the work presented so far, it should be obvious that neither of these error levels is 
immediately relevant in the simulation of the image formation. What is rather more relevant 
is the error incurred on the diffraction orders that are collected by the imaging optics. If the 
spatial frequencies of the error depicted in Figure 12(c) are sufficiently high, this error will be 
discarded by the projection lens. This concept is shown in Figure 13, where the spectra 
(magnitude) of the scattered fields from edge-DDM and from rigorous 3D simulation are 
shown in Figure 13(a) and (b) respectively. The marginal circles indicating the propagating 
plane waves (with \k\<2^/X) and the collected by the imaging system plane waves (with 
1^1 < 2;riV^(l + <x)/m, for NA=0.75. R=4 and a=0.3, are indicated on these plots. The 
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magnitude of the error in the spectrum of the field obtained from edge-DDM is shown in 
Figure 13. Observe that most of the error is concentrated in the higher spatial frequencies for 
\k^\ approaching 0.03 and only small levels of error exist in the collected orders. The 
normalized mean square errors of all propagated and collected orders are 2.2% and 2E-4 
respectively. 

[00681 Similarly, the amplitude of the scattered fields across the observation plane below 
the isolated island obtained with the edge-DDM and with fully rigorous 3D simulation and 
the amplitude of the error are shown in Figure 14. Again, a locally significant error 
amplitude near the comers is evident. However, the energy of this error is concentrated 
mostly at the higher spatial frequencies, for |^,| approaching 0.03 and the error in the 
collected orders is seen to be small, only 0.16% in a NMSE sense. Note that the large DC 
component of the isolated island mask layout was subtracted in order to make meaningful 
comparisons of the edge-DDM with fiiUy rigorous simulations. Otherwise this large DC 
component would artificially show an even better accuracy, since the error level in 
comparison with a much sfronger signal (containing the DC term) would be smaller. 
[0069] Each fiiUy rigorous mask simulation of the patterns of Figure 12 and Figure 14 
that was performed with TEMPEST used approximately l.SGbytes of memory and required 
almost 2 days on a 550Mhz CPU. On the other hand, once all of the required 2D edge 
simulation results are available the scattered field can be obtained with the edge-DDM in less 
than 1 sec. The simulation of the isolated edge-diffractions is an off-line procedure. It only 
needs to be performed once and the edge-diffraction results are recycled as many times as are 
necessary for the reconstruction of the scattered fields of arbitrary layouts, hi the above 
examples, it took approximately 5min for each rigorous, isolated-edge 2D simulation. 
[0070] The edge-DDM is used here to generate the scattered field from a three level 
(0790°/270°) alt. PSM. A portion of the arbifrary target layout is shown in Figure 1 1 before 
and after the phase assignment. Rectangle and edge extraction is performed within the 3^m 
by 4yim (IX) coherence window shown and the scattered field obtained with the edge-DDM 
inside the U^im by 16nm (4X) coherence window is also shown in Figure 1 1 for TE and TM 
illuminations at X=l93nm. Note that a fiiUy rigorous 3D mask simulation of a n^xm by 
16^m (4X) layout would require close to 40Gbytes of memory and would only be feasible on 
a multi-CPU architecture. The scattered field is obtained with the edge-DDM 
implementation in less than Imin. This time includes the polygon and edge extraction and 
the field synthesis from pre-calculated edge-diffractions. Therefore, the results of the edge- 
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DDM cannot be immediately compared with a fully rigorous solution, but the accuracy of the 
edge-DDM within the imaging system's collection ability can reveal the inaccuracies of the 
ideal k-mask, where all transitions are assumed perfectly sharp. This is done in Figure 16, 
where the images (at the 30% intensity level) obtained with the k-mask and the edge-DDM 
for two different undercuts (Onm and 50nm) are compared. The system parameters are 
A,=193nm, NA=0.7, a=0.3 and R=4. The inaccuracies of the k-mask model are revealed at 
various locations in the image as seen in Figure 16. The effect on the image of changing the 
undercut from Onm to 50nm can be quickly and accurately evaluated. Observe the subtle 
differences evident in the intensity contours of the Onm and 50nm masks of Figure 16. 
[0071] Although the edge-DDM presented earlier achieves excellent accuracy and 
tremendous speed-ups as compared to a fully rigorous mask simulation, there is still one 
more parameter to exploit in order to further speed-up the method. There is no need to carry 
all the details of the true field scattering as far as the imaging process is concerned (with 
larger than 1 reduction factors), since the higher spatial frequencies are not used. If a lumped 
parameter model can be devised that closely matches the TTL spectrum of the scattered field 
then the task of an accurate field representation is accomplished. A rect-function, a raised 
cosine function and a Gaussian function can successfully mimic the TTL spectra of scattered 
fields from various size and profile openings. A similar principle can be applied here with 
the purpose of matching the TTL spectra of edgediffracted fields. Assuming a restriction to 
piecewise constant models, the scattered field from an isolated edge under TE or TM 
illumination can be approximated using the multistep fimctions shown in Figure 17. 
Although not shown in Figure 17, the amplitude, phase and size of all steps are optimized to 
achieve the best spectral matching of the collected (TTL) orders with the orders resulting 
from the continuous edge-diffracted field. Note that the key factor for a good specfral match 
is to capture the finite extent of the "off-on" edge-transition. A piecewise constant fimction is 
probably not the best choice to achieve this task, but on the other hand a piecewise constant 
fimction offers a good lumped parameter approach and it lends itself nicely as an extension of 
the unmatched Kirchhoff approach, where all edge-transitions are assumed perfectly sharp. 
Instead of approximating the true edge-diffraction by a single step fimction (geometrical or 
Kirchhoff approximation) a multi-step piecewise constant function is employed. The 
ampUtudes, phases and sizes of all steps are optimized to achieve the best spectral matching 
of the collected orders with the orders resulting from the continuous edge-diffracted field. 
The matched bandwidth approximations to the field diffracted by a 180° edge with 50nm 
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undercut under TE and TM illuminations are shown in (a) and (b) respectively. 
[0072] If the matched bandwidth edge-DDM is applied on the 400nm by 400nm (4X) 
isolated hole of Figure 12 the scattered field across the observation plane can be represented 
in the compact way shown in Figure 18(a). Recall that the edges parallel to the y-axis "see" 
TE polarization and the edges parallel to the x-axis TM polarization. The multi-step matched 
BW approximations to the edge-diffraction shown in Figure 17 for TE and TM illuminations 
intermingle in the x- and y-directions to produce the multi-color mask of Figure 18(a), where 
each color represents a different amplitude and phase region of the mask. The spectrum of 
the error of this approach is depicted in Figure 18(b). Note that compared to the error shown 
in Figure 1 3(c) this error has increased by ahnost four-fold. Nevertheless, the error in the 
TTL spectrum still remains at low enough levels for accurate imaging simulations. A better 
approximating function, other than the piecewise constant step function, that is more capable 
of capturing the true edge-diffraction profile (such as linear or higher order polynomial) 
should result in better spectral matching and consequently more accurate imaging 
simulations. 

[0073] Consider now a domain decomposition strategy appropriate for handling 
wavelength size phase defects. The bright future for aU. PSMs as a major resolution 
enhancement technique (RET) has been plagued primarily by the difficulty to rehably inspect 
them, flag locations that phase defects are present and subsequently repair the defective 
locations. State-of-the-art inspection systems utilize a focused laser beam that can be either 
at-wavelength (up to ?.=248nm), meaning that the wavelength of the inspection beam is the 
same as the wavelength that the alt. PSM was designed for, or emit light of a larger 
wavelength. This beam scans the entire mask and the reflected or transmitted scattered fields 
from the particular mask location under inspection are imaged with a very high resolution 
imaging system that typically has a reduction factor of R=l and NA close to 1. But here is 
the caveat: After the signal from the inspection tool for a particular mask location is available 
it needs to be compared to something in order to determine whether or not a phase-defect (or 
other defect) is present within that location. Such a comparison is in general very difficult, 
because not only does it need to be accurate in order to reliably flag defective locations, but it 
also needs to be rapid enough for the entire reticle to be inspected in a reasonable amount of 



time. 



[0074] The most promising direction that has been adopted by researchers working on the 
problem has been what is known as a die-to-database comparison. By that, it is meant that 
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the local signal is compared to the signal that the same mask location would produce if it 
were defect-free. The problem now shifts to building the database of signals from defect-free 
masks. Note that all possible mask geometries that will be encountered in every inspected 
reticle need to exist in the database. One way to building the database would be to produce a 
(set of) test reticle(s) that includes every single geometry situation that is anticipated to exist 
in all designed layouts to come. This (set of) test reticle(s) is then meticulously inspected 
with a reliable tool (maybe an atomic force microscope - AFM) to guarantee that it is defect- 
free and finally all signals from the different locations are collected and placed in the vault 
(database). Clearly, the immensity of this task has to do not only with the a priori precise 
anticipation of what layout situations to take provision for, not only with manufacturing the 
test reticle(s), characterizing it and building the database, but also with the sheer volume of 
data that such a database would contain. Moreover, a precise diagnosis of the type, size and 
location of a phase defect requires an even larger database of benchmark signals, where every 
possible combination of defect type, size and location in every possible layout configuration 
needs to exist in the database. Performing such a task in a way similar to the one outlined is 
impractical. 

[0075] Simulation can again come to the rescue, at least in principle. If a simulation tool 
that can rapidly and accurately predict the expected signal resulting from a defect-free or 
defective location on the mask is available, then the task of building the database is simplified 
fremendously. Moreover, if such a tool is really fast (and accurate) then a database is not 
even needed. The required benchmark signal from the defect-free layout is generated in-situ, 
while the inspection system gathers measurements. The problem with this solution is that to 
date there exists no simulation tool that possesses both the rapidity and accuracy properties. 
The family of simulation tools that rigorously solve Maxwell's equations around the reticle 
(such as TEMPEST) and subsequently use a vectorial formulation for the image formation to 
calculate the expected optical signal from the inspection system are impractical because of 
large memory and time requirements. On the other hand the accuracy of speedier simulation 
programs that circumvent the solution of Maxwell's equations around the mask is 
unacceptable. 

[0076] Yet another important consideration of an inspection system for aU. PSMs should 
be the following: Phase-defects that are not critically affecting the image that a projection 
printing tool will produce on the wafer surface should be discarded, i.e. not flagged and not 
repaired. But this is also a difficult task, since the only bullet-proof way of achieving this 
would be to expose wafers with the suspect defect present and examine if the printed resist 
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images have intolerable artifacts. Clearly, such a process now involves the close cooperation 
of the inspection system with the projection printing tool. However, an accurate and properly 
calibrated simulation tool can alternatively assess the tendency of the defect to print or cause 
other image artifacts. Consequently, the costly step of exposing wafers can be avoided and a 
decision of repairing or not the defect rests fully upon the simulation result. Again, to date, 
results from no simulation tool, that is fast enough in calculating expected images with phase 
defects present, can be trusted to base repair decisions on. 

[0077] The domain decomposition methods described above are appropriate for solving 
the problem of the rapid and accurate evaluation of the benchmark signal from nondefective 
masks for the inspection system. A domain decomposition method is presented that 
addresses the problem of the defect printability assessment, after a defect is previously found. 
The generation of benchmark signals from defective mask locations can also be based on the 
domam decomposition technique, although more work will be needed to adapt it and test it 
for the high-NA inspection optics with R==l . As will soon be obvious, the speed and 
efficiency of the method are inherited by virtue of the domain decomposition methods 
presented above, although for simplicity edge-DDM is not explicitly employed in the 
foUowmg. In the simulation examples, the illumination wavelength is 193mn and the phase- 
wells are designed accordingly, so as to provide the required phase shifts for that wavelength. 
[00781 Consider the example depicted in Figure 19. In Figure 19(a) the near scattered 
field below a small layout of a 0°/180° alt. PSM is shown when a square 120mn x 120mn 
(4X), 180° post phase defect is present at various locations within the layout. Note that 
although neither the layout nor the phase defect geometry change, a different simulation 
altogether is needed because of the different relative positions of the defect within the layout. 
Each one of the nine 3D simulations in the matrix shown in Figure 19(a) takes approximately 
lOhrs on a 450MHz CPU and utiUze 300Mb of memory. Subsequently, the aerial images for 
all nine situations are calculated for an imaging system with R=4, NA=0.75 and a=0.3. Then, 
the fractional change of CDaeria, at the middle and at the end of the line because of the defect 
is found for all nine defect locations and is shown in Figure 19(b). From the plot in Figure 
19(b) it can be deduced that when the defect is located in the middle of the phase wells or is 
mostly "tucked" underneath the absorbing layer (columns 3 and 1 of Figure 19(a) 
respectively) the CDaeHa, is not critically affected, but locations such as those in column 2 of 
Figure 19(a) cause larger CDaeHa, variations. In the former cases the defects could probably 
be discarded whereas in the latter repair seems necessary. Clearly, the luxury of the lOhr- 
long simulations is only available as a proof of concept and such an approach is not viable for 
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full-chip characterizations. 

[00791 Now, consider the paradigm shown in Figure 20. A domain decomposition 
method is again invoked, where instead of a one-step simulation of the defective layout two 
separate rigorous simulations are performed, one with just the defect-free layout and one with 
just the defect in clear surroundings. The scattered field below the defect is then shifted 
accordingly, such that the defect is effectively placed at the location it appears in the 
defective mask, and the uniform background light is subtracted, so as to retrieve the 
(complex) signal of the local perturbation that the defect causes. This perturbation is added 
to the defect-free layout signal for an approximation to the signal from the rigorous 
simulation of the defective mask. The non-defective mask simulation can be performed 
rapidly using the edge-DDM decomposition methods described above. 
[00801 A sfrategy can now be devised, where a database of the near electromagnetic field 
scattered from isolated phase-defects in clear surroundings is created for all possible potential 
defects that are critical. The volume of such a database can be reduced based on the earlier 
observations that defects with similar footprints (lateral dimensions) and small shape 
perturbations are equivalent from an electromagnetic point-of-view. 
[00811 The suggested decomposition method will be introduced through the following 
example: Suppose a 150mn (IX), 1:1.5 semi-dense, 90°/270° alt. PSM contact mask. The 
layout and the scattered field across the observation plane under normal incidence By (TE) 
polarization are shown in Figure 21. 

[0082] Now, suppose that a 200mn by 200mn (4X) 90° post defect is present in the center 
of the bottom left hole. According to the paradigm of Figure 20, the mask scattering 
simulation can be broken up into two constituent parts. This approach is shown in Figure 22. 
In Figure 22(a) the phase defect is simulated in an isolated configuration at the location that it 
appears in the layout, m Figure 22(b) and (c) the near field across the observation plane 
below the defective alt. PSM is shown, using the domain decomposition method and the 
complete mask respectively. The differences are not discernible from the plots of Figures 
22(b) and (c) and one has to look at the amplitude of the error in Figure 22(d). Again, not 
only is the error level of the decomposition low. but the frequency content is very high, so 
that it will be filtered out by the imaging system. 

[00831 Next, suppose that the same 90° post defect is present in the center of the bottom 
right hole. The important observation here is that no new simulation is necessary for the edge 
decomposition method. The scattering of the isolated phase defect can be recycled from the 
database, i.e., from Figure 22(a). However, the scattered field is now properly shifted in the 
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lateral x-direction and also (a key step) it is propagated in the z-direction. This results in the 
plot of Figure 23(a), where a larger diffraction spreading is evident. In Figures 23(b) and (c) 
the near field across the observation plane below the defective ah. PSM is shown, using the 
domain decomposition method and the complete mask respectively. The amplitude of the 
error is shown in Figure 23(d) and similar observations apply in Figure 22(d). 
[0084] The near fields across cut-lines passing through the center of the bottom two holes 
(90°, 270°) are shown for the two defective masks in Figure 24, comparing the rigorous (r- 
mask) and the decomposition methods (DDM). The non-defective case is also plotted for 
comparison. 

[0085] Next, the images of both rigorous and DDM approaches, for both defective masks 
are compared in Figure 25, (a) for the 90° hole and (b) for the 270° hole. THe optical system 
parameters are X=193mn, R=4, NA=0.75 and a=0.3. The agreement between the r-mask and 
the DDM is good, a normalized mean square error of less than 0.3% is incurred (measured 
only in the area of the defective hole). The k-mask model is insufficient for capturing these 
effects and is only shown for comparison. 

[00861 The reason of the excellent agreement in the predictions of the aerial image can be 
traced back to the plots of Figure 22(d) and Figure 23(d). The high spatial frequency 
variation of the error means that it is mapped at the extremities of the spectrum of 
propagating plane waves and is not collected by the optical system. The normalized spectra 
(magnitude) of the near field errors are shown in Figure 26. Observe that the impUed 
accuracy from these plots should be enough even for R=l and NA~0.8-0.9, which would be 
sufficient for inspection simulations. 

[00871 Thus the invention includes a suitable extension of the domain decomposition 
framework for the simulation of alternating phase shift mask with phase defects. The value 
of this approach is believed to be significant in the problem of rapid assessment of defect 
printability once a defect has been found, or in the die-to-database comparisons of inspection 
systems when deciding about the existence or not of a phase-defect. 
[00881 The invention is an extension of knovm domain decomposition methods that 
introduces tremendous versatility, since only a small number of isolated edge-diffraction 
simulations is shown to contain all the necessary information for the synthesis of the scattered 
field from arbitrary 2D mask layouts and subsequent accurate imaging simulations. The 
limits of this method are reached when the mask features are smaller than a wavelength in 
size and the vertical mask topography is large. Through a systematic process it was 
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determined that features as small as at least 200mii bearing 270° of phase-wells can be 
accurately decomposed via the edge-DDM, at ?.=193nm. The method was tested for simple 
2D layouts, where rigorous mask simulations are possible. Excellent accuracy accompanied 
by speed-up factors of 172,800 (Isec vs. 2days) have been demonstrated. The accuracy of 
the method is attributed to the fact that on one hand the edge-diffraction phenomena are 
modeled rigorously through 2D edge-diffraction simulations and that on the other hand the 
errors incurred during the synthesis (primarily near the comers) are mapped at the exfremities 
of the spectrum of propagating plane waves and do not contribute to the image formation. It 
was emphasized that all necessary rigorous 2D simulations of the diffraction from isolated 
edge profiles is performed off-line and recycled for the diffraction calculation of arbifrary 
layouts. If needed, the accuracy of the edge profile diffraction simulations can be pushed to 
extreme limits, since it is performed only once. This can aid for example in the correct 
simulation of unusual edge profiles whose geometrical details would require excessive 
discretization of the domain. It was also shown that the exact details of the edge-diffraction 
are not necessary for accurate imaging simulations (witii larger than 1 reduction factors) and 
a piecewise constant, multistep edge-transition model that matches the TTL spectrum was 
infroduced. The range of validity of the edge-DDM is expected to be appropriate for rapid 
and accurate evaluation of aerial images whenever speed is critical, as for example in fiiU- 
chip OPC software and die-to-database comparisons in the inspection of alt. PSMs. 
[0089] While the invention has been described with reference to specific embodiments, 
the description is illustrative of the invention and is not to be construed as Umiting the 
invention. Various modifications and applications may occur to tiiose skilled in the art 
without departing from the true scope and spirit of the invention as defined by the appended 
claims. 
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